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Prominent approaches to quantum gravity struggle when it comes to incorporating a positive cos¬ 
mological constant in their models. Using quantization of a complex SL(2,C) Chern-Simons theory 
we include a cosmological constant, of either sign, into a model of quantum gravity. 


The universe is accelerating in its expansion. The sim¬ 
plest consistent explanation for the empirical data is a 
positive cosmological constant A = 10“®^ m“^. A com¬ 
mon expectation is that a quantum theory of gravity will 
shed light on the surprising value of this constant [1- 
3] and on the interplay between quantum fields and the 
bare value of A [4, 5]. A solid evaluation of this proposal 
requires the inclusion of A into quantum gravity. 

Prominent approaches to quantum gravity, such as 
the AdS/CFT framework of string theory, covariant loop 
quantum gravity, and group held theories, struggle to in¬ 
corporate a positive cosmological constant, if they are 
able to incorporate one at all. A few frameworks, like 
causal dynamical triangulations and asymptotic safety, 
do not have this difhculty. 

In this letter we present a novel model of four¬ 
dimensional Lorentzian quantum gravity based on a dis¬ 
cretized path integral over gravitational holonomy vari¬ 
ables where the cosmological constant emerges geometri¬ 
cally as a consequence of our quantization procedure via 
complex Chern-Simons theory. In particular, both signs 
of A are treated on an equal footing in the model. 

Our discretization of the path integral decomposes 
spacetime into a simplicial complex, with each simplex of 
constant curvature A. We choose to work with parallel 
transports along closed paths (holonomies) as they are 
the most natural gravitational observables. They also 
ht nicely with the use of constant curvature simplices, 
whose geometry, as shown below, can be completely en¬ 
coded in a finite number of holonomies. From the Chern- 
Simons perspective, these holonomies arise as the non- 
contractible cycles of a 3-manifold obtained by removing 
a graph from the 3-skeleton of the simplicial complex. 

Before delving into the details of the model, a few 
comments on discreteness, symmetry, and the continuum 
limit. Discreteness is first introduced into the model as a 
regularization tool, as in a lattice gauge theory. However, 
the lattice will be dynamical and fluctuating in a quan¬ 
tum mechanical fashion. We show below that to avoid 
ambiguities in the model the spectrum of the fluctuations 
must have a Planck-scale discreteness and be gapped. 
Hence, the presence and absence of a lattice site is also 


subject to quantum fluctuations. This is the quantum 
‘creation’ and ‘annhilation’ of the grains of spacetime. 

Because of the quantum dynamical nature of these dis¬ 
crete grains, the problem of recovering general relativity 
in the continuum limit is intertwined with a detailed un¬ 
derstanding of their collective semiclassical behavior. 

We do not treat the continuum limit here, however we 
provide evidence for a connection with classical gravity. 
This is achieved by studying the small h asymptotics of 
the quantum amplitude of a single building block. In 
this limit, we recover the exponential of the Einstein- 
Hilbert action (with the appropriate Gibbons-Hawking- 
York boundary term) for a simplex of constant curva¬ 
ture A. This is the curved Regge action of simplicial 
gravity. The bulk action J — 2A) corre¬ 

sponds, on the simplicial side, to a term g/g AU, where 
V is the 4-volume of a Lorentzian and homogeneously 
curved 4-simplex. Intriguingly, both the 4-volume and 
the boundary term are a result of the stationary phase 
evaluation of the Chern-Simons action subjected to ap¬ 
propriate boundary conditions. 

The use of curved, instead of flat, simplices as building 
blocks serves a twofold dynamical purpose. Firstly, the 
interior of each building block is a solution of the dynami¬ 
cal equations of gravity with a cosmological constant and, 
secondly, by gluing a collection of these blocks one can 
build global solutions with the correct continuum sym¬ 
metries [6], that is, those of (A)dS. In this sense, these 
building blocks are maximally adapted to the symmetries 
of the problem and, as such, are promising for studying 
perturbations around a background that includes a cos¬ 
mological curvature—a potentially important technique 
for addressing the smallness of A. 

A further consequence of using these building blocks is 
that their curvature A enters as an infrared cutoff scale. 
Thus, in this model the cosmological constant plays a 
dual role to the Planck area, the latter providing a nat¬ 
ural ultraviolet cutoff. Together these scales make the 
quantum theory of a fixed discretization finite. 

The model presented in this letter has a rich set of re¬ 
lations with other physical and mathematical theories. It 
can be seen as a generalization of the Turaev-Viro [7], and 
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FIG. 1. A planar projection of the tubular neighborhood of 
the Fs graph with the set of longitudinal and transverse paths 
used to calculate the holonomies Gab and Ha{b), respectively. 
For simplices, the Poincare dual is also a combinatorial sim¬ 
plex. 

hence Chern-Simons-Witten [ 8 ], quantization of three- 
dimensional gravity. The physical picture draws heavily 
from loop gravity, especially in its covariant, spinfoam 
version. On the other hand, the mathematical tools come 
largely from the connections that Chern-Simons theory 
has with knot polynomials and the volume conjecture. 
The idea of a relation between Chern-Simons theory and 
four-dimensional loop gravity with a cosmological con¬ 
stant dates back to the early days of these theories [9- 
13]. String theorists have worked extensively on complex 
Chern-Simons theory and we see here potential for ex¬ 
change of ideas between loops and strings. 

In particular, the holomorphic blocks that play a key 
role in defining the present model have an M-theory in¬ 
terpretation, and closely relate to supersymmetric gauge 
theory, as well as surface operators with junctions [14- 
18]. While not holographic in any standard sense our 
model does exhibit a rich interplay between three- and 
four-dimensional topological and gravitational theories 
(for recent work on this interplay in lower dimensions, 
see e.g. [19]). We hope that recent advances in these 
active fields will accelerate progress in quantum gravity, 
for example in treating the smallness and continuum limit 
problems. 

Geometry from holonomies —Parallel transport along 
a Faraday-Wilson loop of the gravitational field results in 
a holonomy that encodes aspects of the curvature of the 
region bounded by the loop. Consider an homogeneously 
curved, geodetic 4-simplex in an (A)dS space, and the set 
of holonomies along loops encircling the simplex’s faces. 
Choose these loops to share a common base point and 
be contained in the 1-skeleton of the simplex. Not all of 
these loops are topologically independent and hence the 
holonomies are subject to a set of algebraic constraints. 
Before displaying these constraints explicitly we move to 
a dual description of the simplex. 

The Poincare dual of a 4-simplex is a graph, we call 
it Pa, that associates a vertex to each tetrahedron of the 
simplex and a graph edge to each pair of neighboring 
tetrahedra. These edges pierce the shared faces of the 
paired tetrahedra and hence are dual to these faces. Pic¬ 


torially one can think of this as concentrating each of the 
faces on distributional defects (i.e. the edges of the dual 
graph) living in a three-sphere topologically equiva¬ 
lent to the boundary of the 4-simplex. With this image in 
mind it is clear that there is an isomomorphism between 
the fundamental group of the simplicial 1 -skeleton and 
that of the graph-complement manifold M 3 = S^XNfT^), 
where A^(r 5 ) is the open tubular neighborhood of the 
graph. 

Therefore the loops used to calculate parallel trans¬ 
ports are in 1 -to-l correspondence with a set of equiva¬ 
lence classes of loops in the graph complement M 3 . As¬ 
sociating the parallel transport holonomies to the corre¬ 
sponding equivalence class of loops we obtain a represen¬ 
tation of a flat connection in M 3 . The boundary of M 3 , 
or equivalently of the closure of A^(r 5 ), is a genus six 2 - 
surface, call it Se. Any loop of M 3 can be deformed to 
live on this boundary. An (overcomplete) basis of non- 
contractible loops on Ee is given by a set of 4 x 5 loops 
transverse to the tubes and 10 lines running longitudi¬ 
nally along the tubes, see Fig. 1. Labeling the vertices 
of Fs with a = 1,..., 5, the longitudinal holonomy from 
a to 6 is Gba = G~b, and the transverse holonomy based 
at a and encircling the tube ab is Hb{a). 

In this basis the topological constraints mentioned 
above break up into three classes. They are: 

closures flt,Mb(a) = 1 ; ( 1 ) 

parallel transports GbaHb{a)Gab = ( 2 ) 

, ,, f j GacGcbGba = ^, 

bulk contractions < , , (3) 

|G'34G'42G'23 = 

The closures are a consequence of contractibility around 
each vertex. The parallel transport conditions are due 
to contractibility around each tube. The last class of 
equations come from contractibility in the bulk of M 3 and 
are expressed in terms of a set of six independent cycles. 
The first five are (abc) G {123,125,134,235, 345}, see Fig. 
1, while the last is explicit. Eq. (3) is a consequence of 
the single essential crossing in the planar projection. 

Returning to these holonomies as gravitational observ¬ 
ables on the 4-simplex we must account for the fact 
that each tetrahedron of the simplex defines a three- 
dimensional frame. This is why the four holonomies 
{Hb{a)}b associated to tetrahedron a should stabilize the 
direction transverse to this frame (when appropriately 
parallel transported [20-22]). We take the transverse di¬ 
rection to be timelike and therefore these holonomies are 
part of an SO(3) subgroup of the Lorentz group. When 
viewed as a condition on the fiat connection of M 3 we 
call this requirement the ‘geometricity constraint’. 

We see that corresponding to each 4-simplex geometry 
in (A)dS there is a flat connection on M 3 satisfying the 
geometricity constraints. In a generalization of a classic 
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theorem due to Minkowski [23, 24], we have shown that 
the converse also holds: 

Theorem 1. There is a bijection between a dense set of 
Lorentz flat connections on satisfying the geometric- 
ity constraints and non-degenerate convex constant cur¬ 
vature 4-simplex geometries, up to parity. [21, 22] 
Because this is a bijection, a flat connection can be used 
to reconstruct either a positively or a negatively curved 
geometry, but only one of the two. For a precise state¬ 
ment of the non-degeneracy requirements, see [ 20 - 22 ]. 

Sketch of the proof: The SO(3) little group associated 
to each vertex identifies the timelike transverse direction 
to each tetrahedron. Forming the appropriately paral¬ 
lel transported scalar products among these directions 
we can build the matrix of cosines of the hyperdihedral 
boost angles 0ab- This Gram matrix fully encodes the 
geometry of a unique 4-simplex. One finally checks the 
consistency of the reconstructed geometry with the other 
data contained in the holonomies, such as the areas ao& 
of the triangular faces. The determinant of the Gram 
matrix encodes the sign of the curvature. 

Chern-Simons quantization —Having translated the 
space of four-dimensional simplicial geometries into the 
space of flat connections on a three-manifold, it becomes 
natural to quantize this space using Ghern-Simons (CS) 
theory. Indeed, CS theory is the only three-dimensional 
gauge theory with flat connections as solutions to the 
equations of motion. To connect with the geometrical 
reconstruction above we still miss one ingredient, the ge- 
ometricity constraint. We impose this constraint as a 
boundary condition on the CS path integral on M 3 . 

Both for mathematical convenience and with the in¬ 
tention of eventually coupling this theory to fermions, 
we lift the Lorentz group to SL(2,C). With I[A\ = 

/ms \^AdA -\- §41^] the complex action is, [25, 26], 

S[A,A] = '^I[A]+^-I[A], (4) 

where A is an SL(2,C) connection in the fundamental 
representation, k is the complex level (or inverse cou¬ 
pling), and the overbars indicate complex conjugation. 
Efforts to define precisely the SL(2,C) theory have, in 
part, relied on analytic continuation of SU(2) CS the¬ 
ory [27], see also [28, 29]. For our present purposes, the 
formal path-integral definition is sufficient. 

With quantization in mind we seek a set of canonically 
conjugate coordinates in the classical theory. A standard 
and convenient choice are complex Fenchel-Nielsen (FN) 
coordinates x and y, which can be understood as follows 
[22, 30, 31]. We restrict attention to the holomorphic sec¬ 
tor, with the other one obtained by complex conjugation. 
The Atiyah-Bott symplectic structure, H, of CS theory 
sets transverse and longitudinal parts of the connection 
as conjugate. A symmetry reduction of the gauge trades 
the connection variables for our transverse and longitu¬ 
dinal holonomies [32], which also turn out to be conju¬ 


gate. More precisely, if we introduce the eigenvalues of 
the transverse and longitudinal holonomies 

a^ab = and = (5) 

respectively, then {uabiVcd} = Sab,cd- To complete the 
variables introduce a trinion (pants) decomposition of 
each 4-valent vertex and the associated {ua,Va}a- 

The variables {ua, Va}a are purely imaginary as the ge- 
ometricity constraint imposes that they are in an SU(2) 
subgroup of SL(2,C). The reconstruction theorem im¬ 
plies that Uab oc iA&ab, where a.ab is the area of the corre¬ 
sponding triangle in the 4-simplex, and hence Uab is also 
purely imaginary. Finally, ‘^Vab is complex and its real 
part is oc Oab, which is the hyperdihedral boost angle 
between tetrahedra a and b. All of the logarithmic co¬ 
ordinates u and V have the standard branch ambiguity, 
leaving their values defined only up to an element of 27rfZ 
and this will play a role briefly. 

With the classical coordinates specified we turn to con¬ 
structing the physical wave function of a quantum grain 
of geometry. We do this in the u-representation and 
for notational convenience introduce u = (ua,Uab)- The 
wave function is constructed using the path integral 

Z{u,u)= f VAVAc^^^AA]^ (g) 

J u,u 

where the subscripts indicate that the integration is per¬ 
formed at fixed u. This wavefunction satisfies a set of 
Schrodinger-type equations, A(i, y; k)Z{u, u) = 0 and its 
anti-holomorphic analog, which are quantum mechanical 
implementations of the topological constraints of Eq. (3) . 
As functions of the FN coordinates, the A are known as 
A-polynomials, and are a specific quantization of the cel¬ 
ebrated topological invariant of 3-manifolds [33-35] . The 
wave function Z(u,u) factorizes, [14], 

Z{u, u) = T,ana,Z°‘{u)Z°‘{u). (7) 

Here, Z°‘{u) and Z“(it) are the holomorphic and anti- 
holomorphic 3d-blocks, and a labels different branches 
of the solutions to the A-polynomial equations. The re¬ 
construction demonstrates that there are two branches, 
a = 1,2, corresponding to the two parities of the 4- 
simplex. In terms of FN coordinates these solutions differ 
only in the sign of the real part of ^Vab, that is, in the 
sign of the boost Qab- Requiring suppression of other 
possible branches can be used to fix the coefficients 
and is under study. 

Note that the holomorphic block Z'^{u) is defined 
through an analytic continuation of the SU(2) CS theory, 
hence the name ‘holomorphic’ [14, 27, 35, 36]. While the 
Darboux coordinates u and v are useful for the WKB 
analysis below, it turns out that non-perturbatively the 
holomorphic block is a meromorphic function of x and 
thus must be periodic in u. This is connected with the 
finiteness of our model discussed in the conclusion. 
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FIG. 2. A sketch of the Lagrangian manifolds £a and and 
their intersections in un-space. 


WKB analysis —To understand the semiclassical, 
small h, limit of these wavefunctions we use a geometri¬ 
cal approach to the WKB approximation. A WKB wave- 
function is of the form ■ The phase function 

la is a solution of the Hamilton-Jacobi equation and is 
most easily constructed using the geometry of Lagrangian 
submanifolds of the phase space [37]. Here we focus on 
the relative phase between the two branches of Eq. (7). 
We comment on the overall phase below. 

To calculate the phase difference I 12 = I 2 — h, we 
have to choose a loop £ in phase space connecting the 
two critical points and integrate the Liouville form along 
it, I 12 = Jgi9. For the result to be a solution of the 
Hamilton-Jacobi equation, the loop is chosen to lie along 
two Lagrangian submanifolds. The first one, Ca., is im¬ 
plicitly defined by the equations of motion A = 0. The 
second, fixes the representation of the wavefunction 
and is simply the plane u = const. These two manifolds 
intersect precisely at the stationary phase points, where 
a pair of parity related classical geometries is picked out, 
see Fig. 2. 

For calculating the integral of the Liouville form the 
choice of path on a Lagrangian manifold is locally irrele¬ 
vant due to the defining property H|£ = 0. One remain¬ 
ing subtlety is the shift of the phase as you integrate 
through a caustic of £a, where the geometry is degener¬ 
ate, this is the well-known Maslov index 77 [38, 39]. We 
defer study of this index, which may vanish [25, 26, 34], 
to future work. 

An effective technique to calculate Ii 2 {u) is to consider 
small deformations of the loop. To this end, we consider 
two nearby solutions differing by a small amount 5u. The 
variation 5Ii2 splits into contributions along Lu+Su, 
and Ca- The first two contributions vanish because 1 ? = 
udu. Thus, 

Kk 

5Ii2 = E0at,(5aab (8) 

and is calculated along the manifold of geometrical solu¬ 
tions Ca- The 5ua contributions cancel between the two 
branches. Hence this has exactly the form of the Schlafli 


variation [40, 41]. The Schlafli identity then allows us to 
conclude that 

I 12 = Y 2 -Tii ~ ^ 12 ) I (9) 

where C 12 is a geometry independent integration con¬ 
stant, possibly including 77 , that we drop in what follows. 
The geometric origin of the Schlafli identity guarantees 
that the sign of the A in the action agrees with that of the 
reconstructed 4-simplex. Restricting to the real contour 
of functional integration, where A and A are complex 
conjugated variables, we obtain 




cos 


AS(fc) 

Vlirh 


AH 


( 10 ) 


up to a global phase, which simplifies in the case of 4- 
simplices in the bulk of the triangulation. Identifying 
the coefficient with we recognize the 

curved Regge action calculated on our semiclassical grain 
of spacetime. This result is analogous to the asymptotics 
of the Turaev-Viro model for three-dimensional quantum 
gravity. 

While discussing the branches of Ca we ignored the 
lift ambiguity of the logarithmic FN coordinates. When 
properly accounted for, this ambiguity enters the argu¬ 
ment of Eq. (10) and adds a term 

( 11 ) 

bn 

where the Nab are arbitrary integers. Therefore to have 
an unambiguous state we require this expression to be an 
integer multiple of 27r. This forces a quantization of the 
areas in units of SnjGh, where 7 = 9(fc)/3fi(fc). 

This concludes the WKB analysis of a single 4-simplex. 
Just as a simplicial complex can be built from gluing in¬ 
dividual 4-simplices we can glue graph complement mani¬ 
folds using surgery at the graph vertices; the geometricity 
conditions must be imposed at all surgery sites. At fixed 
areas the WKB analysis goes through simplex by simplex 
and gives a quantum simplicial complex weighted by its 
global curved Regge action. A full understanding of the 
physics on the whole complex and its continuum limit is 
involved. For example, what happens when you sum over 
the areas (the flatness problem) and the role of the local 
fluctuation of a simplex’s parity (the cosine problem) are 
recognized difficulties under investigation. 

Crucially, although the geometrical reconstruction pro¬ 
ceeds simplex-by-simplex, the cosmological constant and 
its sign must be semiclassically consistent across the com¬ 
plex. A WKB solution exists only if the curvature of all 
the simplices is the same. Note that this does not mean 
that the simplicial manifold cannot support curvature 
defects as in Regge calculus. This consistency is a conse¬ 
quence of the fact that the sign of the curvature is already 
determined at the level of the tetrahedra in the boundary 
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of the simplex and therefore propagates to neighboring 
simplices. 

Relation with the EPRL model —A natural way to im¬ 
pose boundary conditions on the connections in M 3 = 
S^\T5 is to insert an appropriate Fs Wilson graph op¬ 
erator in 5^. A Wilson graph operator satisfying our 
desiderata can be built out of the EPRL intertwiners at 
the core of spinfoam quantum gravity [42, 43]. These 
desiderata are the geometricity constraints and their con¬ 
sistent implementation with quantized areas. In the spin- 
foam context the geometricity constraints are referred to 
as ‘simplicity constraints’ and are used to turn topolog¬ 
ical BF theory into general relativity by imposing that 
the R-field is the square of the tetrad field, B = e Ae. 

The area spectrum follows from the quantization of 
the classical variables via SU(2) and SL(2,C) representa¬ 
tion theory. This procedure associates the spin quantum 
numbers of SU(2) to geometrical areas. The spins of 
SU(2) are half-integers, not integers, and this introduces 
an area-dependent sign in Eq. (10), already present in 
the A = 0 spinfoam asymptotics [44] and in the Ponzano- 
Regge model [45] . The spectral spacing is determined by 
the Planck area, SttGH, multiplied by the dimensionless 
Barbero-Immirzi parameter, and the latter can be iden¬ 
tified with our 7 . 

Early work [46] on loop gravity attempted to quantize 
the theory with 7 = 1 but ran into difficulties [47] . Ana¬ 
lytically continued CS theory allows for a continuation of 
both 3fi(fc) and 3(fc) to arbitrary complex numbers [27]. 
This freedom provides an opportunity for spinfoam mod¬ 
els, since it might be used to consistently continue the 
real Barbero-Immirzi parameter back to the imaginary 
unit. This continuation is related to the original self¬ 
dual Ashtekar variables, with multiple advantages such 
as a more intimate relation to covariant geometries, a 
simplification of the dynamics, and compelling black-hole 
entropy counting [48, 49]. 

The present work gives a CS theory generalization 
to four dimensions of the quantum group deformation 
used in three-dimensional quantum gravity models 
to implement a non-vanishing cosmological constant. 
Although for three-dimensional models finiteness on 
a fixed simplicial complex is limited to the positive 
cosmological constant case [7], in the model presented 
here finiteness extends to both signs of A. In fact, 
finiteness is a consequence of the presence of a finite 
number of states, see [50, 51], which in turn follows from 
the periodicity and discreteness of u, Eqs. (5) and (11). 
This illustrates one of the, we hope, many advantages of 
including a cosmological constant in quantum gravity. 
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